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Abstract 

. . . Multi-qudit systems are studied in tomographic probability representations of quantum 

^ ■ qudit states. Results of calculations for Bell-type numbers within the framework of classical 
Q ! probability theory and in quantum tomography are compared. Violation of Bell-type inequali- 
CN I ties are shown explicitly using the method of averaging in tomographic picture of quantum states. 
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p^; 1 Introduction 

C ' The Bell-CHSH inequalities pLj, L2j and their violation are considered as specific properties of 
^ ■ non-local quantum correlations. On the other hand, recently [3], [4] it became clear that the 
O^i Bell inequalities can be associated with properties of some probability distributions in classical 
I probability theory f5]. Prom this point of view violation of Bell and Bell-type inequalities [6], [7], 
CN ■ tSj, tSj can also be studied in the framework of standard probability theory. The aim of this work 
1^ . is to extend the results of our previous work [10] and establish Bell-type inequalities for random 
ly-^ I systems with specific joint-probability distributions. We show that there exist families of joint- 
'sj" ■ probability distributions depending on extra parameters which violate Bell-type inequalities. 

. We apply tomographic probability representation [TT], [T2| of quantum qudit states and show 
^ I that the joint-probability distributions (tomograms) describing the multiqudit states can violate 
' Bell-type inequalities. The systems of two qubits, three qubits, two qutrits and qubit-qutrits are 
O ■ considered explicitly. The paper is organized as follows. In Section 2 we present the tomographic 
^ I probability representation approach to studying multipartite systems. We discuss the origin of 
Bell-type inequalities and their relation to quantum non-local correlations. The violation of Bell- 
type inequalities is linked to the general problem of separability and entanglement of quantum 
5^ , states. We illustrate the theory on the well known example of two-qubit system and point out 
several cases of separable and entangled states based on the analysis of Bell-type inequality for 
the system and its upper bounds. In Section 3 we study the qubit-qutrit system in both classical 
and quantum tomographic probability representations of multipartite states. We establish both 
classical and quantum upper bounds for Bell-type inequalities in qubit-qutrit system. Several 
examples of parameters violating classical Bell-type inequalities are presented and their domains 
are calculated. Finally, in Sections 4 and 5 two-qutrit and three-qubit states are studied in terms 
of establishing the upper bounds of possible Bell-type inequalities analogous to Tsirelson bound 
in two-qubit case. 

A set of non-negative real numbers normalized to unity can be considered as some proba- 
bility distribution P{n) with a finite number of outcomes. One can then rearrange this set 
to receive yet another probability distribution P{m). One can say that P{m) is induced by 
P{n). Moreover, some of these induced probability distributions can be interpreted as marginal 
distributions for the initial P{n) which in its own turn can be regarded as a new joint-probability 
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distribution Wij. out tnen having a joint probability distribution one can associate it witn tne 
probability distribution of observables pertaining to measurements in corresponding subsystems 
and construct various functionals in the form of average values of these observables. However, 
it is clear that the fact of existence of such marginal and joint-probability distributions strongly 
depends on the properties of the initial set or table of probabilities. Thus, if they indeed ex- 
ist, the functionals constructed as described above must by necessity provide some inequalities 
similar to Bell-CHSH or analogous constraints on joint-probability induced functions. It is in 
this context that one of the aims of our work is to show that and how Bell-type inequalities are 
associated with probability distributions or stochastic matrices deduced from the latter. The 
violation of Bell-type inequalities by Quantum mechanics is due to the specific properties of 
quantum tomograms inducing joint-probability distributions with special properties. 

Our work follows conceptually the idea of possible generalization of Bell-type inequalities 
over quantum systems with arbitrary number of qudits similar to |13| . However, the method 
of tomographic description of quantum states [H] allows one to constitute and examine higher 
dimensional systems based only on the properties of some specially constructed probability 
vectors and stochastic matrices. The variety of possible Bell-type inequalities in a system made 
up of s subsystems of di, d2, ■ ■ ■ dg dimensions accordingly, is conditioned by the choice of some 
special sign matrix denoted by our convention as Idiisid2^...^ds- We develop a phenomenological 
method for constructing this matrix by considering the order of signs entering the expressions 
of both classical and quantum-wise averaged observables. In principle, the approach could allow 
to unite many concurring Bell-type inequalities within the framework of a standard probability 
theory. Through our work on the example of two-qutrit and three-qubit systems families of 
Bell-type inequalities are discovered which are not violated by quantum tomograms of entangled 
states. This results echo back to the well-known "openness" of the problem of establishing the 
widest constraints on correlations by local realistic theories [E], |Tn|. It then turns out that Bell- 
type inequalities and their upper bounds B^ax of multi-qudit states can be written in a general 
form by introducing a single correlation matrix C generalizing the special matrix Idx®d2»...®ds- 
This allows to view Bmax as some general functions defined on the elements Cmi,m2,...,ms of the 
matrix C. This result is quite valuable to us as one can then construct multi-moment correlation 
functions in search of stronger Bell-type inequality violations. Throughout our work we see that 
if the Bell-type inequalities of an entangled multiqudit state are violated they tend to violate 
stronger than for two-qubit system [T7], although the factual search of such states is in practice 
much more complicated, too. A general survey and discussion of Bell-type inequalities was 
presented in [T8| . 

2 Observables and Correlations 

Let's examine the general approach of calculating averages and correlations in multipartite sys- 
tems within the classical tomographic representation of states. For simplicity consider composite 
system of two qubits or two-level subsystems. In the past few years several new Bell-type inequal- 
ities were discovered ^ to pinpoint non-local variable theories violated by quantum mechanics. 
Our method has the additional advantage of tightening or loosing the Bell-type inequalities by 
the appropriate selection of a special matrix C which is inspired by recent findings of tight 
Bell-type inequalities [T3] . Consider qubit as a quantum-mechanical spin 1/2 system with only 
two internal degrees of freedom. A quantum-mechanical observable C as a function of local 
observables which in this case are the spin projections m can be represented by a 2 x 2 matrix 
in the standard spin 1/2 basis |m),m = ±1/2. Its quantum-mechanical average can be written 
in terms of classical probability theory: 

{C) = Cmm'W^,m' = XjC V) ("^1^^ V') = Tr (c^mL^^] = {C\mL^'^) (1) 

m,m' m,m' 

- (-) 

Here we have introduced the 2x2 stochastic matrix M^^ built up from tomographic prob- 
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abilities ■w[+-^,u ) and ■w[—2,iJ ) oi qubit up (m = +2) down" [m = —■^) states, u is 

the unitary rotation matrix entering the definition of the quantum tomogram |19| . [20] . From 
the other hand, given the density matrix p of the quantum state one can write for the average 
of C the well known quantum-statistical formula: 



{C) = Tr{pC) = {p\C) 



(2) 



A substantial difference between the two ways of determining the same average value (C) is 
that in the case of (2) density matrix contains probabilities entering the final expression only via 
its diagonal elements whereas the matrix is wholly built of probabilities. It is true that in 
case of a solitary qubit the difference is purely mathematical since there are only two probabilities 
w{+^) and w{—^) entering the final expression for the average value. However, in case of complex 
systems built up of several qubits the first method (1) (i.e. the method based on determining 
the average of observables via stochastic matrices and probability vectors associated with the 
state p) has a clear advantage over the method based on density matrix description of states as 
it allows the experiment at or to view all the quantities involved as physical observables i.e. the 
matrices C, quantum tomograms and M^,. Moreover, if one wishes to measure average values of 
observables in a rotated coordinate systems in (2) one needs to conduct SU (2) transformation 
for both p and C i.e. write 

p^U^pU, C ^f{U)Cf{U), 
where T{U) is a unitary transform acting on the original C matrix. 

In quantum tomographic representation (1) one only needs to indicate different directions in 
space entering the tomograms as column- vectors in matrix (see below the matrix (6)). This 
allows to construct Bell-type inequalities and observe their violation in principle at an arbitrary 
number of directions. These are directions along which measurements can be carried out. Thus, 
to discover new Bell-type inequalities and upper bounds of multiqudit states one simply needs 
to set the correlation matrix C and construct the stochastic matrix Myj. Let's illustrate the 
above said on the example of a system made up of only two qubits. 

In two-qubit system one registers four possibilities for the composite state Wj^ , w |_, 

w which can be arranged into a stochastic matrix: 



w 



W ^ W \- W ^ W ^ 

w w w w 



(3) 



The correlations between dichotomic variables mi and 1712 inside the composite system are given 
then by formulas 



(mim2) 



mi ,m2 



mim2 ■ Wmima = Tr [a2^2Mw ) = {cr2m\M, 



(4) 



where we have introduced the scalar product of operators in the standard spin 1/2 basis \mi) (gi 
\m2). The matrix a2^2 consists of unity elements taken with the same signs with which each of 
the probabilities Wij enter the expression for correlations. It reads: 



C2iS)2 



(5) 



Since w++ + W-i +w |- + w = 1, for the considered stochastic matrix (3) one will always get 

the inequality \{mim2)\ < 1 which is a simple Cauchy-Schwartz inequality that the correlation 
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cannot exceed i m absolute value, in [6) tne tour columns nave tlie interpretation oi discrete 
probability distributions. In order to predict non-local correlations, we must equip the theory 
with physically meaningful parameters which can be measured experimentally. Most commonly 
these are the Euler angles setting directions for two qubits to be oriented in space. Thus, 
varying the probability column- vectors entering (3) along some four fixed directions x, y, z, t 
and assigning them different letters we arrive at a "realistic" stochastic matrix: 



x,y,z,t 



Xl 
X3 



yi 
y2 
2/3 



Z3 



h 

t2 



1 - Xl - X2 - X3 1 - yi - y2 - 2/3 I - Zl - Z2- Z3 l-ti-t2-t3 



(6) 



Each of variables Xj, yj, Zk, ti, i, j, /c, / = 1, 2, 3 in (6) can assume values between and 1. The 
columns of this stochastic matrix represent four possible probability distributions of the two- 
qubit system. These probability column-vectors correspond to two distinct space orientation of 
each qubit in the composite system. In case of quantum probability distributions the variables 
Xi, Uj, Zk, ti become functions of 16 Euler angles and thus render into continuous probability 
distribution^. 

The results (3) — (6) allow generalization of our method in calculating correlations inside mul- 
tipartite systems. The quantum-mechanical average of an observable C describing an expected 
correlation between local observables mi and 1712 (with matrix elements Cmim2 = ("^i|C'|"^2)) is 
given by the formula: 



E 



'^mi m,2 ,m'j^ rTij m.2 rrtj 



J2 {m[m2\C^ 



mi ,1712 ,m^,m2 



\mi ,m2 



1 " (--) 



Tr (^C^m!) 2 



{C\M, 



where the measurements of each of local observables m, m' can yield up to {2d -|- 1) x {2d' + 1) 



outcomes, d 



m. 



max 1 d 



From these definition of correlations eligible both in the classical local (hidden) variable 
theories and the quantum mechanics it follows that an arbitrary linear "correlation form" is the 
familiar average value: 



5(x,y,z,t) =rrCMx,y,z,t 



(7) 



where the matrix Mx,y,z,t is the stochastic matrix (6) associated with the two-qubit system. 
The word "linear" refers to our way of calculation in standard probability theory as expressed 
in (1). The form of (6) is the most general expression for a stochastic matrix associated with 
two-qubit states both in classical and quantum cases. The building blocks for the theory are the 
matrices ddi_'^d2®...®ds and Idi®d2®...®ds as concrete realizations of a general correlation matrix 
Cdi®d2®...<sds violating Bell-type inequalities. The mechanism of constructing Ofii^d2®...®ds was 
described in [TO]. The matrix Idx(^d2®...®ds is a simple induction of the latter by transposition 
and sign inversion operations and was also defined there. The multivariate function B (x, y, z, t) 
reminds one the function that results in any of various Bell-CHSH inequalities and may be 
in that respect further named Bell-CHSH typ^ form. It was established in [10], [12] that in 
classical local variable theories the Bell-CHSH form is a harmonious function of its arguments. 
This formally results in Bell-type inequalities in classical probability theories. It is then natural 
to call the maxima of the function B (x, y, z, t) Bell-type numbers. In this paper we call the 



^In reality there are six Euler angles attached to each tomographic probability distribution entering the resulting 
stochastic matrix as its columns. However, the azimuthal angles V'n, n = 1, 2, 3, 4 don't enter the final expressions. Thus 
the total number of angular parameters necessary to determine stochastic matrix associated with a given composite 
quantum state p reduces from 24 to 16. 

^For simplicity we will be sometimes referring to the form (7) as simply Bell-CHSH form. 
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greatest m value among all-possible r5ell-type numbers tlie upper-bound oi tlie multi-qudit state 
of study. In composite systems the Bell-type numbers (and upper-bounds) are a unique but not 
complete measure of the strength of both local classical and non-local quantum correlations. 

If the composite state is separable it contains no non-local correlations and the Bell-CHSH 
form is limited to some maximal value. In some cases the Bell-type numbers may even frequently 
vanish. Indeed, if we choose the matrix C equal to a constructed from all-possible combinations 
of dichotomic variables mi and m2 we notice that: 



B (x, y, Z, t) = rrd-2(xj2Mx,y,z,t = ^Xi - ^Vj - ^^k + ^ti 







However, if the system is entangled the joint-probabilities entering the stochastic matrix Mx,y,z,t 
don't factorize and the non-local quantum correlations may overthrow the classical upper bounds 
or maximal correlation strengths. 

It is interesting to note, that by the special choice of matrix C the Bell-CHSH form reduces 
its dependence only to four parameters. Indeed, the direct calculation shows 



± 7 

B (x, y, z, t) = rrPMx,y,z,t = 2 (xi + yi + zi 



-ti-l) 



where the matrix P differs from an earlier matrix l2^2 from | ilQj with the change of sign in the 
last column. Thus in the case of matrix P the composite state is described by the same number 
of parameters as in the case of separable two-qubit state. Had we taken the matrix a2<^2 in form 
of (5) we would have received all Bell-type numbers null just as the matrix associated with the 
joint-probability distribution would have shown. The Bell-type inequalities and upper bounds 
are received by examining the properties of the B-form in (7). The minimums and maximums 
of the Bell-CHSH form are achieved only at the 2^ = 16 boundaries of probability distributions 
Xj, yj, Zfc, ti in full agreement with an earlier similar approach described in [12]. All sixteen 
Bell-CHSH type inequalities are restrained in the classical domain of < 2 i.e. the upper 
bound in case of two qubits is 2. In Section 7.1 of Appendix the results of analytic and numeric 
computations for two qubit system are presented and here we would like to stop on detailed 
description of the analysis of the quantum case of entangled system of two qubits. 

The tomogram of state described by a 2 x 2 density matrix p is presented by the probability 
vector 



w (p, U] = Diag (u^pU 



kllPpil + UUU21P2I + U21U11P12 + \u2l\^P22 



\ui2?Pll + U12U22P2I + U22U12P12 + \U22?P22 



Entering rotation operator matrix from SU (2) in the general form, we arrive at 



w{p,e,(p) 



pn cos 



sm 



9Re {pi2( 



-I- P22 sin^ 



2 

6 
pii cos^ - + sin 6Re (/9i2e~*'^) + P22 cos^ - 



where and 9 are the Euler angles entering the expression for rotation operator for spin -. We 

see that the qubit tomogram does not contain Euler angle explicitly. 

In \\^ we found that the upper bound of Bell-type number within the framework of classical 
probability theory applied to the two-bit system is 2. For the maximally entangled state 



1/2 








1/2 


























1/2 








1/2 



(8) 
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we observe violation oi Jjell-Uilbil inequality tor two-qubit system at a local point and near tne 
upper bound as summarized in Table 1. Following the methodology of tomographic probability 
theory we choose two directions 1,2 and l',2' for each of the qubits and predict quantum 
non-local correlations between their respective dichotomic variables. In our theory these non- 
local correlations rise naturally as an effect of stochastic nature of the matrices constructed of 
tomograms measured at various combinations of chosen qubit directions. 



Angles 


Directions 


Directions 




1 2 


r 2' 


1 2 


r 2' 




6.2800 6.2800 


6.2800 0.0416 


6.2639 1.8010 


0.0190 0.0200 


e 


1.5834 2.3705 


2.3705 0.7962 


1.5834 3.1410 


2.3620 0.7910 




B = 2.4148 


B = 2.8284 



Table 1: Violation of Bell-CHSH inequality for qubit-qubit system near the Tsirelson bound. 



Taking the values of Euler angles from Table 1 we are able to indicate the stochastic matrices 
resulting in the violation of Bell-CHSH inequality. 



Mb-. 



:2.4148 



0.4264 0.4263 

0.0735 0.0736 

0.0735 0.0736 

0.4264 0.4263 



0.4999 0.2490 

0.2509 

0.2509 

0.4999 0.2490 



Mb-. 



:2.8284 



0.4280 


0.4255 


0.4278 





0742 


0.0720 


0.0745 


0.0722 





4258 


0.0720 


0.0745 


0.0722 





4258 


0.4280 


0.4255 


0.4278 





0742 



In Section 7.1 of the Appendix the calculation of the explicit expression for 
B2»2iP2(s,2,0l,^Pl,62,^P2,63,V3,^4,Vi\C = 12^2) are given. 

The second column of directions in Table 1 corresponds to the well known case of so called 
Tsirelson bound [2l]- It manifests to the maximum entanglement in the qubit-qubit system 
in terms of the strongest quantum non-local correlations between the pair of qubits; a result 
which is beyond the scope of any classical local variable theory. The value of the upper bound 
is symmetric as it follows from the analytic expression of the Bell-CHSH form (see Section 7.1 
of Appendix). Numeric computation allows to obtain the lower bound as well: 



91 = 1.9244 

92 = 0.5657 



1.9244 = 1.7351 993 = 1.7351 
0.5657 ip2 = 0.7780 lpa^ = 0.7780 



B = -2^/2 

as expected the result is symmetric as it should be for any general linear correlation strength 
restriction. 

In Fig. 1 the results of calculations of Bell-type numbers as functions of the angle 61 are 
shown. The picture graphically demonstrates Tsirelson's theorem [21]. In Section 7.1 of the 
Appendix the surface of Bell-CHSH form is drawn near the Tsirelson bound of two correlating 
angles 02 and ^4. 

In the next several sections we study multiqudit states of higher dimensions {d > 2) and 
calculate Bell-type numbers in classical probability theory. We then calculate Bell-type numbers 
and upper bounds of Bell-CHSH forms in the quantum tomographic theory and compare results. 
Our aim, as clarified in the Introduction is to point out examples of multiqudit states where 
Bell-type inequalities exist with various upper bounds analogous to Tsirelson bound in the lower 
dimension of d = 2. With the implementation of quantum tomographic methods we see that 
there exist whole domains of parameters violating these Bell-type inequalities. 

Our approach allows to easily generalize the results obtained for the maximally entangled 
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1.4 



0.5 



1.5 



2.5 



3.5 



Figure 1: The Bell-CHSH form i?2®2(^ii V^ii ^2, '^2, ^s, '^3, ^i, ¥^4) in the vicinity of Tsirelson bound as 
a function of 9i. 



state (8) over well-known Werner states for two qubits [22] ■ Indeed, for the state 
1 +P 



p 
2 





1 — p 
4 







4 




2 




1 +P 



1 ~ p 
4^ + 1 



1 
— 

' 1 

10 0- 



The Bell-CHSH form assumes the expression: 

B^2{Gl,^l,92,^2,03,ip3,04,ipi,p) = 

(OiV fd^V /OiY fe^V 

= 2p-4pCOs( yj -4pCOS( yj -|-4pCOslYj COSIyj -|- 
+4pCOs(yl COslyl -F4pCOs(yl COS ( y I -4pCOs(yl COS ( y I + 



+'ip COS 



9i\ f93\ 
— cos — 
2 J \2j 



+4p C0S( y I COSI y 
+Ap C0S( y I COsI y 
-4p C0S( y I COsI y 
-Ap COS( y I COSI y 



'^™(^) C0S(933)-F 

sin^y^ sin^y^ cos(v9i) cos(9J4) -h 
sin( yj sinfy j cos((/?2) cos(993)- 



sm| y j sinf y ) cos((/?2) C0S(934) 



sm 
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-4pcos^yJ '^°®\^tJ ^™VtJ ^™VTJ sinCy^i) sin((^4) - 
-Apcos^^^ '^°^(y) ^™(y) ^^"(y) sin(<^3) + 

+4pcos^y^ *^°^(y) ®^^(y) ^^^(y) ^^^'^'^2) sin(</?4). 

Figure 2 shows the dependence of upper bounds of Bell-inequahty for two qubits on the parameter 
p that counts for the entanglement of the system. As it can be drawn from the graphic, the 
upper bound grows hnearly in the interval of entanglement ^ < p < 1 reaching its maximum at 
the value of p = 1 . 



3 




-S.4 -0.2 0.2 0.4 0.6 0.8 1 



P 



Figure 2: The upper bounds of Werner states as function of parameter p. 



The two-qubit system reviewed in this chapter provides an insight into studying higher dimen- 
sion multiqudit states at least with respect to the amount of information that the Bell-CHSH 
form of a multipartite system stores in it. First of all the tomographic probability representa- 
tion of multipartite states equips us with a universal tool both in classical and quantum realms. 
Procedure- wise the program of studying the system's Bell-type inequalities and upper bounds can 
be formulated the following way. Each subsystem contributes to the resulting joint-probability 
distribution of the composite system via a stochastic matrix. Information of correlations between 
various dynamic variables in the system (both in quantum and classical cases) are contained in 
the stochastic matrix of the composite state which can be either simply separable, separable or 
entangled (non-separable). This matrix enters a bilinear form named Bell-CHSH form through 
the average value of the product of correlating observables. We calculate all local (Bell-type num- 
bers) and global maxima (the upper bound) of the Bell-CHSH bilinear form. We establish first 
the Bell-type numbers in classical probability theory for multi-level systems. Then we conduct 
the calculation of Bell-type numbers through quantum tomograms and compare results. This 
recipe allows one to detect possible Bell-type inequalities, classify them by their corresponding 
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upper bounds and predict violation oi tnese inequalities m quantum mecnanics by observable 
angles rising in the measurement process of the quantum system. Moreover, given an entangled 
multiqudit state our method allows to find explicitly the domains in the parameter space where 
Bell-type inequalities are preserved and the small quantum domain where they are maximally 
violated. We will now realize this program for several interesting multiqudit states. 



3 Bell-CHSH type inequalities in qubit-qutrit system 

The correlations inside the qubit-qutrit system with two chosen directions ni and 112 are given 
by the expectation value 



{17111712, 



ni, n2 



mi =±1,7712 =±1,0 



w {mi, 1712, 111,112) rnim2 



(9) 



Just as in qubit systems the orientations of qudits in space are brought in by the necessity of 
measurements. In principle, one can choose any number of possible polarizations for each of the 
qudits in the composite system. However, since the signature of the correlation matrix of the 
bit-trit system can be presented in the form of a 6 x 6 matrix 



a(2® 3) 



h- ' 


1 


1 


1 


1 


1 




















-1 


-1 


-1 


-1 


-1 


-1 


1 

1 L 


-1 


-1 


-1 


-1 


-1 




















1 


1 


1 


1 


1 


1 



(10) 



we prefer ordering three directions of qutrit and two directions for qubit. This allows one to 
conduct maximal analysis of the resulting B-form. The expression for special sign matrix (10) 

is obtained from formula (9) by assuming all probabilities equal to -. The special signature 
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matrix 1(2 ® 3) is received by transposing a{2 3) and inverting the sign of its last ro^MJ. 

To calculate the upper bounds of Bell-type inequalities i.e. obtain all possible Bell-type num- 
bers, we need to consider 2^ = 256 possible stochastic matrices related to 5-possible directions, 
two per two-level system and three per the three-level system. However, not all of these vertices 
are physically feasible since the half of these vertices contain values of parameters Xi, yj, z^, ti 
{i,j, k,l = 1, 2) prohibited by the requirement of positiveness of stochastic matrices i.e. 



1 - yi - y2 > 0, 



1 - zi - Z2 > 0, 



i-u 



t2 > 0. 



so the number of eligible vertices is reduced to 2^ x 3^ = 108. The two-level subsystem and the 
three-level subsystem are described by two stochastic matrices respectively 



Xi 
1 — Xl 



X2 
1 - X2 



y,z,t 



1 



yi 
y2 

yi - y2 



1 



Zl 

Zl - 



Z2 1 



tl 
t2 

tl-t2 



:ii) 



where each of the variables Xi, yj, Zk, U {i,j,k,l = 1,2) assume values between and 1. The 
stochastic matrix of the composite system is Mx,y,z,t = Afx "X" -^y,z,t and the calculation for the 

first 27 vertices for B-form B (x,y,z,t) = Tr (C'Mx,y,z,t) gives all-possible Bell-type numbers 



^Generally speaking, various techniques of receiving various matrices / from the generic a matrix are irrelevant 
to the problem of entanglement and subsequent violation of Bell-CHSH type inequalities. As long as the difference 
between the quantum and classical Bell-type numbers (upper bounds) are calculated via the same matrix C the 
developed approach stays perfectly legal. 
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tor the qubit-qutrit system m the general torm as linear tunctions oi the elements oi the matrix 
C . The list of Bell-type numbers in case of these vertices is given in the appendix. 
The corresponding Bell-CHSH type numbers for a particular matrix 



1 





-1 


-1 





1 


1 





-1 


-1 





1 


1 





-1 


-1 





1 


1 





-1 


-1 





1 


1 







-1 





1 


-1 





1 


1 





-1 



(12) 



are calculated to be 

B = -6, 0, 6, 6, 4, 4, 2, 2, 0, -4, 0, 2, 6, 8, 12, 6, 0, 2, 2, 4, 4, 6, 0, 2, 0, 0, 6 

The values ±6, ±8 and ±12 correspond to non-physical vertices (see formula (1) in Section 7.2 
of Appendix for clarification) and should be discarded. Thus, we obtain 

TDclassical a 

^max ^ 

Consider now the quantum case. Let's examine the qubit-qutrit entangled state 



P2ig)3 — 



1/2 














1/2 










































































1/2 














1/2 



The classical bound \B\ = 4 can be reached on many points in the ten-dimensional phase 
space spanned by ten Euler angles, four for the two directions of qubit and six angles for the 
three directions of qutrit. One of the examples that can be instantly checked is: 



Anglos 


Directions 


Directions 




1 2 


V 2' 3' 


1 2 


1' 2' 3' 










0.03 3.19 





e 





TT TT TT 


0.03 0.22 


3.08 3.09 1.65 




D = -[ 


B = 1.0312 



Table 2: Violation of Bell-CHSH inequality for qubit-qutrit system 



The corresponding stochastic matrix easily calculates at Euler angles entering Table 2 and 
gives: 
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1 
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4 


4 


4 
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The violating parameters for qubit-qutrit system were found through searching very small 
quantum domains. One example of such violation is presented by the second column of Table 
2. An example of a stronger violation of Bell-type inequality for qubit-qutrit system is given in 
the Table 6 of Section 7.2 of the Appendix. 
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The stochastic matrix Mx,y,z,r,s,t of the two three-level classical system reads: 

Mx,y,z,r,s,t = M>c,y,z Afr,s,t 

where each of the matrices entering the final expression has the form of (11). The explicit 
expressions for Bell-CHSH form and Bell-type numbers as functions of an arbitrary correlation 
matrix Cs^s can be found in Section 7.3 of the Appendix. 

There are 2^^ = 4096 possible vertices and corresponding Bell-type numbers in this case. 
Again, because of the requirement of the positiveness of probabilities entering the stochastic 
matrices, this number is reduced to 3^ = 729 physically eligible vertices. 

The upper bound corresponding to the particular case of the correlation matrix I303 is found 
to be (see formula (3) in Section 7.3 of the Appendix) 

ryclassical y 



5 Bell-CHSH type inequalities in three-qubit system 

We consider now the system of three qubits. The classical probability theory approach treats 
the system as an eight-level system. The first 14 Bell-type numbers as functions of arbitrary 
matrix elements (720202 are presented in Section 7.4 of the Appendix. Number of vertices i.e. 
number of possible Bell-type numbers in this case is 2® = 64. Using the standard /2(gi2(gi2 matrix 
we obtain the following Bell-type numbers in classical tomographic probability theory: 

B = 0, ±2, ±4, ±6 

The upper bound of Bell-type numbers is therefore 

■Ticlassical ry 

^max 

In interesting discussion of three qubit system correlation functions was indicated in [25. In 
Section 7.4 of Appendix explicit expressions of classical and quantum correlation functions are 
presented in form of B-forms. These expressions change dramatically depending on the values of 
matrix elements of the sign matrix C. By running a numeric search for local Bell-type inequality 
violations for several particular cases of C one needs to exceed the value B = 6. The quantum 
violation of Bell-type inequalities here are very subtle due to large amount of possible inequalities 
not all of which can be violated. 



6 Conclusion 

To conclude, we formulate our main results. We studied several multiqudit systems in terms of 
their possible Bell-type numbers restricting correlations in multipartite states. Within the clas- 
sical framework of probability theory we used simple tomographic approach to describe classical 
states as probability distribution vectors. In the quantum case the tomographic approach was 
sophisticated by entering quantum tomograms continuously depending on many Euler angles 
and strongly violating classical Bell-type inequalities. Bell-type numbers were constructed as 
maxima of some B-forms arising naturally in our tomographic description of states as average 
values of non-local correlation between subsystems. Our results are not complete as in case of 
three-qubit and two-qutrit systems the choice of the special matrix C becomes quite difficult and 
the violation of Bell-type inequalities becomes possible for many matrices. It is possible that 
there is no one universal matrix I for theThis matter needs a further analysis which we plan to 
present in our future works. The results of Bell-type numbers received in the article is presented 
in the Appendix with detailed description of numeric and analytic aspects of calculations. It is 
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notable tnat the values oi derivations oi quantum mecnanics trom classical local variable tlieories 
based on probability description of both classical and quantum states depend on the same ma- 
trix C (or /) i.e. there are hundreds of local variable theories violated by quantum mechanics. 
Partial comparison of these theories are expressed in terms of Bell-numbers in classical local 
variable theories and in quantum mechanics. Some of our results include: 

Bbit-bit = 2, Bbit-trit = 4 

Bqubit-qubit = 2.8283, Bqubit-qutrit = 4.0612 

In our next work we plan to present the results for Bell-type numbers in two-qutrit and and 
three-qubit systems for several entangled states including the Werner family of states i.e. one- 
parametric families of Bell-type inequalities. Because the Bell-type numbers only partially char- 
acterize the phenomenon of entanglement in our next works we are going to introduce also various 
information entropies for multiqudit states hoping to establish (or prove otherwise) possible link 
between Bell-type numbers and entropies. In particular we are interested in domains where indi- 
cations of entanglement are manifested both by maximal violation of Bell-type inequalities and 
relationships between entropies of subsystems with entropy of the composite system. We would 
also like to aim at obtaining numeric and analytic results for a system of qudits within classi- 
cal local variable theories to investigate further the origin of special matrix C strongly limiting 
possible Bell-type numbers. We intend to study in future works the multi-moment correlation 
measures, in comparison to those that consider only 2nd moment (pairwise or quadratic studied 
in this paper) dependence. One of the interesting possibilities here is constructing matrix C 
through correlations like 

Ccx (mf _..m^^) 

As it is well known, using multi-moment correlation functions it is possible to get a measure 
for more general dependence of Bell-type numbers on the initially given data (multi-variable 
states, correlation matrices, orientations in space etc.). We expect to increase the numeric 
"gap" between classical and the quantum Bell-type numbers. 

In the present work the numeric calculations were carried out in MATLAB. 

7 Appendix 

In Appendix we have collected useful analytic and numeric results which demonstrate the calcu- 
lation schemes for classical and quantum tomographic approaches. For convenience. Appendix 
is broken into subsections one per each of the systems studied. 

7.1 Two-qubit system 

The classical tomographic probability theory|T2j applied to the system of two-bits provides us 
with Bell-type numbers as functions of elements of an arbitrary correlation matrix C. The 
Bell-CHSH form for two qubits with a fixed matrix 12^2 has the form 

B2«,2ixi,X2,yi,y2) = 2 + 4xi2/i -h 4xi7/2 - 4xi - 4yi - 4x2^2 + 4x22/i. 

The results for 16 Bell-type numbers in case of the familiar matrix /2g)2 are gathered in Table 
3. The combinations of vectors x and y represent the vertex of the polytope in the probability 
space of qubits. Bell-type numbers are the values of Bell-CHSH form taken at these vertices. 
The maximum value among these numbers, namely i? = 2 is the upper bound of Bell-type 
inequality for two-qubit system in classical probability theory. 
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-2 


01 


10 


2 


10 


10 


-2 


11 


10 


2 


00 


11 


-2 


01 


11 


-2 


10 


11 


2 


11 


11 


2 



Table 3: The Bell-type numbers of qubit-qubit system for the matrix /2(g)2 



The matrix P reduces the dependence of B (x, y, z, t) on 12 parameters in non-separable 
two-qubit state into four parameters as in separable state. It has the form: 



P = 



1 


-1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 


-1 


-1 




1 


1 


1 



To calculate Bell-type numbers in quantum mechanics, we utilize two unitary matrices Ui, U2 
for the first qubit and two unitary matrices Vi, V2 for the second qubit. 
The Bell-CHSH form then calculates: 



= 2 + 4 COS (I) ' COS (I) % 4 COS (I) ' COS (I) % 

"2 \ / "3 \ A I '^^ \ I ^4 \ . / "1 \ A I "3 



+4 COS — COS — — 4 COS — COS — — 4 cos — — 4 cos , 

2; ^2; \2j \2j \2j V2 

-4cOS^y^ COS^y^ ^™(^) ^^^('^l) C0S(<^3)- 

—4 cos -— cos -— sm -— sm -— cos[(pi) cos{(p4) — 



^2 J \2 J \2 J \2 ^ 

-4cOs(^y^ COs(^y^ ^™(y) ^™(^) COs((y33) 
-|-4cOS^y^ ^™(y) C0S((^4) + 

-f4cos^y^ cos^y^ ^'^'^{^^ ^™(t) ^^^("^i^ sin(v73) 
-F4cos^y^ cos^y^ ^™(y) ^^('^i) sm{ipi) + 

-F4cos(^y^ cos(^y^ ^™(y) ^™(^) ^^^'^^j sin((^3) - 
-4 cos^^yj cos^y^ ^^^(y) ^^^(y) ^^'^'^2) sin((/J4) 

In Fig. 3 the Bell-CHSH form -B202 (^i , V'l , ^2, ¥'2, ^3, <P3j ^4, </'4) as function of 62 and ^4 was 

plotted. The values of the rest of the Euler angles were taken very close to those which give 
Tsirelson bound (see Table 1). The white area on the plot indicates the small quantum domain 
where the Bell-CHSH form exceeds the classical restraint of iJ = 2. The global maxima of the 
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^4 



^2 



Figure 3: The Bell-CHSH form -B2(gi2(6'ij ¥^1, 6*2, 932, ^3, <^3j ^4, ^4) as a function of 6*2 and 6*4 

function is tiie Tsirelson bound. The Bell-CHSH form in case of a general matrix C reads: 

B- — + — + — + — + — + — + — + — - 
"2 2 2 2 2 2 2 2 

Oil cos Ci2C0S^^^ CiiCOS^^^ Ci3COS^^j Ci2C0S^^^ CmCOS^^^ 

2 2 2 2 ^ 2 2 ^ 

2 ..o .,0 ..o ..o .0 

Ci3 COS ( ^ j Ci4 cos ( 

2 2 

C22COS^^^ C23COs( 



2 2 2 2 2 2 

OS^^j CuCOS^^j C2lCOS^^^ C22COS^^j C23COS^^j C2lCOS^^^ C24COS^^j 

~2 2 2 ^ 2 2 2 2 

2 2 2 2 2 2 2 

Os(%) C23COs(|) C24COs(f) C3lCOs(f) C3lCOs(|) C32COs(f) C32COs(f) 

^ 2 2 2 2 2 2 

C33COs(|)' C33COs(|)' C;5iCOs(%)' c,iCOs(|)' C4lC0s(^)' C42COs(%)' C4lCOs(f)' 



2 2 2 2 2 

2 



C43COS^^^ C42COS^^^ C44COS^^^ C43COS^^^ C44COS^^^ 

2 2 2 ^ 2 2 

CllCOSIyl COslyl - Cl2 COS ( y 1 COS I y 1 - C13 COS I y 1 COS I y 1 + C14 COS I y 1 COS I y 1 + 

C2lCOSIyj COslyj - C22 COS ( y 1 COS I y 1 - C23 COS I y 1 COS ( y 1 + C24 COS ( y 1 COS I y 1 + 

C3lCOSIyl COSIyl - C32 COS I y 1 COS I y 1 - C33 COS I y 1 COS I y 1 + C34 COS I y 1 COS I y 1 + 

C4I cos I y I cos I y I - C42 COS I y I COS I y I - C43 COS I y I COS I y I + C44 COS I y I COS I y I - 

Cll COs(^y^ ^°^(^) 003(993) + 

C12 COS^yj ^°^(y) ^^"^(y) ^^^('^l) C0S((^3) + 

14 



UJ '°'UJ ''""UJ ''""UJ 

03 



COS 



COS 



sm 



sm 



sm 



sm 



Ci3 COS 
Ci4 COS 
C21 COS 
C22 COS 
C23 COS 
C24 COS 

C31 COS I — I COS I — I smi — I sm 



y j C0S(V71) 



COS 



sm 



— cos — sm — sm 



cos 



sm 



cos 



COS 



sm 



sm 



— cos 



y 1 sm 



C43 cos — COS — sm — sm 



C44 cos — cos — sm — sm 



cii cos 

Ci2 cos 
Ci3 COS 

Ci4 cos 



COS 



COS 



COS 



COS 



C21 COS — COS — sm — sm 



C22 COS — COS — sm — sm 



sm 



sm 



sm 



sm 



C23 COS 
C24 COS 
C31 COS 
C32 COS 
C33 COS 
C34 COS 



COS 



cos 



COS 



sm 



sm 



sm 



sm 



sm 



sm 



— cos — sm — sm — 



sm 



C32 cos 

C33 cos 

C34 COS 
C41 COS 

C42 COS I — - I COS I — I smI — - I sm 



03 



sm 



sm 



— sm 



sm 



sm 



sm 



sm 



sm 



sm 



sm 



sm 



cos(v?i) 
cos((^i) 
cos(<^i) 
cos((^i) 
cos ((^2) 

^ C0S(^2) 

COs((^2) 

COS ((^2) 

C0S(^2) 
y) COS ((^2) 
COs((^2) 

cos(<^2) 

y) sin((^i) 

y) sin((^i) 

y) sin(^i) 

y^ sin((^i) 

y) sin((^i) 

y) sin(^i) 

y) sin((^i) 

y) sin((^i) 

y) sin((^2) 

y ) sin(<^2) 



sm 



W2 

2 

02 

— I cost ^ I sm| — I sm 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



cos 



sm 



sin 



sin 



sm 



sin 



sin 



sm 



sin 



sin 



sm 



y) sin((^: 
yj sin(^2) 



>2) sin 



sin 



(<^3)- 

ifs)- 

{'P4) + 
(¥'4) + 

{^3) + 

{'P3) + 

i'Ps)- 

{'P3)- 
{^4) + 
{'P4) + 

{'P4) + 

i'Ps)- 

i'P3) + 
i'P3) + 
(<^4)- 
ifi)- 
i'P4) + 
i'P4) + 
i'P3)- 

i'Ps)- 

{'P3) + 
{'P3) + 



C41 cos|^yJ ^°^(^yj ^^"^(^yj '^^'^vYJ sin(992) sin(v94) 



C42 COS 



C43 COS 



C44 COS 



COS 



COS 



sin 



sm 



cos 



sm 



sm| — \ sin(992) sin(v94) - 



sm| — J sm((/?2) sin(934) + 



sm| y J sin(932) 8111(994) 



For illustration of the situation when Bell-type inequality is preserved by the two-qubit 
system we present the following values of Bell-CHSH forms gathered in the Table 4. The table 
shows values of the Bell-CHSH form corresponding to two possible orientations for each of the 
qubits denoted 1,2 and l',2' accordingly. 



Anglco 


Directions 


Directions 




1 2 


V 2' 


1 2 


1' 2' 






TT 


71 71 


71 


e 


2 « 




fi 4 


4 4 




B = 1 


B ^ 1.83 



Table 4: Bell-CHSH inequality validation for qubit-qubit system 

Numeric computation also allows to recover the stochastic matrices corresponding to these 
values of Bell-CHSH form. For the value of i? = 1 the stochastic matrix reads: 

1/4 1/4 + 1/8^2 1/2 1/4 1/8^ 

1/4 l/4-l/8\/2 l/4-l/8\/2 

1/4 1/4-1/8^2 l/4-l/8\/2 

1/4 l/4-M/8\/2 1/2 1/4-^1/8^/2 

For an increased but still "classically suppressed" value of B ^ 1.83 its stochastic matrix reads: 



Mb= 



Mb. 



;1.83 



5/16 + l/16\/6 l/4-M/16\/6 3/8-M/16V2 3/8 

3/16-l/16\/6 1/4-1/16^/6 1/8-1/16a/2 1/8 

3/16-1/16^/6 1/4-1/16^/6 1/8-1/16^/2 1/8 

5/16-^1/16^/6 l/4-M/16^/6 3/8-M/16^/2 3/8 



These states are separable and they satisfy the Bell-CHSH inequality |20| . 



7.2 Qubit-qutrit system 

Some Bell-type numbers as functions of arbitrary correlation matrix 6*203 stie organized in the 
Table 5. These numbers are the values of the corresponding Bell-CHSH form at given vertices 
constructed by employing permutations of probability vectors x, y, z,t. 
The Bell-CHSH form for qubit-qutrit system is the function: 

i?2®3(x, y,z,t) = 4 + 4xizi - 2x2*2 + 4x12/1 + 2xiy2 + 4x2^i + 4x2yi + 2x2Z2+ 
+2xiZ2 + 2xit2 + 4:Xiti - 6x1 - 2x2 - 4:yi - 4zi - 2y2 - 2 Z2 + 2 X2y2 - 4 ^2*1 

The quantum Bell-type numbers are calculated by exploiting two unitary matrices from 
SU{2) for the qubit and three unitary matrices from SU{2>) for the qutrit. This produces 
six quantum probability distributions entering the final expression for the stochastic matrix 
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V Lex 
00000000 
00000001 
01000000 
01100000 
01010101 
10000000 
10101010 
11000000 
11100000 



Cl6 + C26 + 
Cl6 + C35 + 
C16 + C26 + 
Ci4 + C26 + 
Cl5 + C35 + 
Cl3 + C33 + 
C21 + Cii + 
Cl3 + C33 + 
Cll + C33 + 



C56 + C46 + C36 + C66 
C26 + C56 + C46 + C65 
C53 + C43 + C36 + C63 
C53 + C4I + C36 + C63 
C25 + C52 + C42 + Cq2 
C23 + C56 + C46 + C66 
C31 + C54 + C44 + C64 
C23 + C53 + C43 + C63 
C23 + C53 + C41 + C63 



Table 5: The Bell-type numbers of qubit-qutrit system 



associated with the entangled state P2i^3- This way the first Euler angles violating the classical 
local reality inequalities(Bell-type inequalities) were discovered. The Bell-CHSH function in 
quantum mechanics is: 

52®3(^l, ^2, ^3, O4, 65) = COs(^i) COS(03) + COS(01) COs(04) + COs(02) COs(03) + 
COs(^l) 008(^5) + 003(^2) COS(04) — 003(^2) COS(05) 

Table 6 summarizes our current results for maximal violation of Bell-type inequalities in 
qubit-qutrit system. The upper bound for Bell-type inequality for the qubit-qutrit system was 
obtained by choosing two directions 1 and 2 for qubit and three directions l',2',3' for the qutrit. 



Anglos 


Directions 




1 


2 


r 


2' 


3' 




0.0004 


3.1450 











e 


0.2368 


0.2538 


3.1471 


3.1292 


1.6024 




^203 = 4.0612 



Table 6: A stronger violation of Bell-CHSH inequality for qubit-qutrit system 



7.3 Two-qutrit system 

The stochastic matrix Mx,y,z,r,s,t associated with the composite system of two trits is taken in 
the form of (10). The elements of special sign matrix as^s are received by arranging the signs 
of probabilities entering the average value 

(mim2)ni,n2 = X] 777,2, ni,n2)mim2 

TOi,m2=±l,0 

into columns of the corresponding sign matrix: 
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-1 
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-1 



















































































-1 


-1 
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ine cell-UMbM lorm m tms case reads: 



-S3®3(x, y, z,r, s,t) = 7 - 6xi - 6yi - 2zi - 3x2 - Sy2 - Z2 - Gn - 6si - 2ti - 3r2 - Ss2 - ^2+ 
+4zisi — 421*1 + '2ziS2 — Z2t2 — 2^2*1 + Z2S2 + 2z2Si + Z2r2 + 2z2ri + +4xiri + 2xir2+ 
+4xisi + 2xiS2 + 4x1*1 + 2x1*2 + 2x2ri + X2r2 + 2x2Si + X2S2 + 2x2*i + ^2*2 + 4yiri+ 
+2yir2 + 4yisi + 2yiS2 + 4yi*i + 27/1*2 + 2y2n + y2?^2 + 2y2Si 
+1)282 + 2 2/2*1 + 2/2*2 + 4 ziri + 2 zir2 — 2 zi*2 

Some of Bell-type numbers for qutrit-qutrit system are gathered in Table 7. These numbers 
were obtained by evaluating the Bell-CHSH form i?3ig,3(x, y, z, r, s. t) at all-possible vertices 
constituted by permutations of x, y, z, r, s, t probability vectors. Using the bilinearity of the Bell- 
CHSH form one can say with certainty that to obtain the upper bound of Bell-type inequality 
one needs to simply compare these numbers for a given matrix 6*303. 



Vertex 


Bell-type numbers 


000000000000 
000000000001 
100000000000 
010101010101 
101010101010 


C49 + C59 + Cgg + C79 + Cgg + C99 + Ci9 + C29 + C39 
C49 + C59 + C68 + C79 + C89 + C98 + Ci9 + C29 + C38 

C49 + C59 + Cgg + C79 + Cgg + C99 + Ci3 + C23 + C33 
C45 + C55 + C65 + C75 + C85 + C95 + Ci5 + C25 + C35 

C51 + cei + C71 + C81 + C91 + cii + C21 + C31 + C41 



Table 7: The Bell-type numbers of qutrit-qutrit system 



This way the several numeric values of Bell-type numbers corresponding to the particular 
case of matrix 6*303 = -^3(g>3 can be obtained: 

B = -9, -3, 0, 6, 9, 15, 18, 6, 0, -12, -18, -30, -36, 0, -6, 0, 3, 9, 12, 18, 12, 0, -6, -18, -24, -12, -6, -9 

The values ±12, ±15, ±18, ±24, ±30, ±36 are non-physical as they correspond to fictitious 
vertices (see page 9) and should therefore be discarded. The corresponding quantum Bell-type 
numbers are received by using six matrices from SU{3) three per each of the correlating qutrits. 
The resulting 9x9 stochastic matrices enter the expression for Bell-type numbers along with 
the special sign matrix /303. These expressions were evaluated at various values of Euler angles 
in the 12-dimensional phase space. The Bell-CHSH form corresponding to sign matrix /s^s has 
the form: 

2cos(gi-g4) ^ 2 cos{0i - 9n) ^ 2cos(g2-^4) ^ 2 cos(gi - gg) , 

-03®3(,Pl, P2, ^3, ^4, P5, f^ej — 2 ' 3 ' 3 ' 3 

2 COS (02 - ^5) 2 008(^3 ~ 94) 2 COS (02 - de) 2 cos(g3 - O5) 2 cos(g3 - de) 
^ 3 ^ 3 ^ 3 ^ 3 3 

It is notable that this B-form does not depend on angles ipi,i = 1, 2, 3, 4, 5, 6 for the chosen form 
of the matrix 1303 and predicts a maximum value B^ppgr = 6 < 7, i.e. the Bell-type inequality 
is preserved for the particular choice of matrix C^^s. 



7.4 Three-qubit system 



We start with following stochastic matrices built-up from probability distributions corresponding 
to each of the 2-level subsystems in the three-bit system: 



Xi X2 
1 — Xi 1 — X2 



yi 
1 - yi 



y2 

1-2/2 



Zl 



Z2 
1-Z2 



The resulting 8x8 stochastic matrix of the composite system Mx,y,z = (X" My iX" has 
the following elements entering its columns Xi2/i2:i, Xi2/i(l — zi), xi{l — yi)zi, Xi{l — 2/i)(l — zi), 
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- xi)yizi, (L - xi)yi[L - zi), [L - xi)(L - yijzi, [1 - xi)(L - yi)[L - zi). iJie second and 
other columns have the same structure obtained with cychc permutation of x, y, z. 
The Bell-CHSH form in the framework of classical theory of probability reads: 

-B2(g)2®2(x, y, z) = -6 + 8x2yiZi + 8xiy2Z2 + 8x2y2Zl - 8X2y2Z2 - 8xiZi + 8x2yiZ2- 

—8yiZ2 — 8yizi + 8x12/2^1 — 8x2^:1 — 8x122 — 8x2yi — 8x11/2+ 
+8x1 + 4x2 + 8yi + 4y2 + 8zi + 4z2 — 8yiXi + 8xiyi^;i + 8xiyiZ2 — 8y2Zi 

The calculation for the first 14 vertices via arbitrary matrix (72^202 yields to following Bell-type 
numbers gathered in Table 8. 



Vertex 


Bell-type numbers 


X 


y 


z 




00 


00 


00 


C18 + C28 + C38 + C48 + C58 + Cqs + C78 + C88 


10 


00 


00 


CU + C21 + C34 + C44 + C58 + C68 + C78 + C88 


11 


00 


00 


Ci4 + C24 + C34 + C44 + C54 + C64 + C74 + C84 


11 


10 


00 


C12 + C22 + C34 + C44 + C52 + C62 + C74 + C84 


11 


11 


00 


C12 + C22 + C32 + C42 + C52 + Cq2 + C72 + C82 


11 


11 


10 


Cll + C22 + C31 + C42 + C51 + Cq2 + C71 + C82 


11 


11 


11 


Cll + C21 + C31 + C41 + C51 + C61 + C71 + C8I 


01 


01 


01 


C18 + C27 + C36 + C45 + C54 + C63 + C72 + C8I 


10 


10 


10 


Cll + C22 + C33 + C44 + C55 + C66 + C77 + C88 


00 


00 


01 


C18 + C27 + C38 + C47 + C58 + C67 + C78 + C87 


00 


00 


11 


Cl7 + C27 + C37 + C47 + C57 + C67 + C77 + C87 


00 


01 


11 


Cl7 + C27 + C35 + C45 + C57 + C67 + C75 + C85 


00 


11 


11 


Cl5 + C25 + C35 + C45 + C55 + C65 + C75 + C85 


01 


11 


11 


Cl5 + C25 + C35 + C45 + C51 + C6I + C71 + C8I 



Table 8: The Bell-type numbers of three-qubit system 
The selection of correlation matrix 



C, 



2(g)2(g)2 



2®2(gi2 



-1 -1 


1 -1 


1 


1 -] 


-1 -1 


1 -1 


1 


1 -] 


-1 -1 


1 -1 


1 


1 -] 


-1 -1 


1 -1 


1 


1 -] 


-1 -1 


1 -1 


1 


1 -] 


-1 -1 


1 -1 


1 


1 -] 


-1 -1 


1 -1 


1 


1 -1 


1 1 


-1 1 


-1 


-1 1 



-1 

as well as choices with sign inversions determine the upper bound of the system by comparing 
all Bell-type numbers. The upper bound is 6. 

The quantum correlations are taken into consideration via the simplest entangled state 
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ine cell-UMbM torm tor the tiiree-qubit system reads: 



-B2(gi2(gi2(^'l, ^3) ^3; ^4, ^6) ^3, ^4, V'Sj V'e) — 

8 C0S(-V32 f6) COS|^y I COS ( y I COS ( y I Sill ( y I Sill ( y I Sill ^ 

, . ^A f^A f^A ■ f^A ■ f^A ■ (^^ 

COS(^-(y9l -ip-i- (/J6j COS( y I COS I y I COS ( y I Sill ( y ) Sill ( — ) Sill ( — 
COS( — — 974 — j COS( — ) COS( — ) COS( y ) Sill ( — ) smj — ) smj — 
COS(-(/?2 - V^S - V^Sj COS( — ) COS( y I COS I y I Sill I y I Sill I y I Sill I y 



-8 cos(-(/?i - 9J4 - fo) ^'^^{^^ 



,'4 \ /6'6 \ . / Ci \ . / (74 \ . ^ (76 

COS — COS — sm — sill — sill 



2 ) \2 J \2 J \2 J \2 

cos[—ip2 — (fs — (fe) cos — cos — cos — sm — sm — sm 



\2 J \2 J \2 J \2 J \2 J \2 

2 \ / ^4 \ / ^5 \ . / ^2 \ . / ^4 \ . / ^5 



cos{-ip2 - 9^4 - V'sj cos( y^ ^^^\2 ) ^^^\2 ) ^^\'2 ) ^^\'2 ) '^"\ 2 

« ( \ f^A f^A f^A ■ f^A ■ f^A ■ f^A 

-8 C0S(-991 -<f3- f5) C0S( y I COSI y I COS ( y 1 Sill I y I Sill ( y I Sill ( y I . 

We hope to find the deviation from the classical bound for Bell-type number due to quantum 
correlations between three-qubits in our next article. The finding of these deviation is compli- 
cated by the necessity of choosing several scenarios of pair interactions inside the system. These 
interactions in form of remarkable correlations between dynamic variables are absorbed in a 
single matrix C with several physical constraints ^0\. The three-qubit B-form obtained above 
in principle permits violation of Bell-type inequality. In that perspective in our approach the 
sign matrix 120202 is much "better chosen" than 130,3 for the two-qutrit system. 
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